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In quantum interferometry, it is vital to control and utilize nonlinear interactions for achieving 
high- precision measurements. Attribute to their long coherent time and high controllability, ul- 
tracold atoms including Bose condensed atoms have been widely used for implementing quantum 
interferometry. Here, we review the recent progresses in theoretical studies of quantum interferome- 
try with Bose condensed atoms. In particular, we focus on the nonlinear phenomena induced by the 
atom-atom interaction and how to control and utilize these nonlinear phenomena. Under the mean- 
field description, due to the atom-atom interaction, matter-wave solitons appear in the interference 
patterns, and macroscopic quantum self-trapping exists in the Bose-Josephson junctions. Under the 
many-body description, the atom-atom interaction can generate non-classical entanglement, which 
may be utilized to achieve high-precision measurements beyond the standard quantum limit. 
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It is well known that an ambitious goal of quantum 
physics is to control and exploit quantum coherence and 
entanglement. Linear superposition of multiple quantum 
states is a perfect description for quantum coherence. 
Macroscopic quantum coherence (MQC), the quantum 
coherence among an ensemble of particles, may lead to 
novel phenomena never existed in single-particle systems 
and new applications in quantum interferometry. The 
systems of Bose-Einstein condensates (BECs) 1-4], such 
as, superfluids, superconductors and Bose-Einstein con- 
densed atoms, are typical matters of MQC. 

Under ultralow temperature close to absolute zero, as 
the thermal motions are almost frozen, the quantum fea- 
tures become significant and robust. In the last 20 years, 
there appear dramatic progresses in our ability to con- 
trol and manipulate atoms. Now, atomic quantum gases 
and atomic BECs can be easily prepared by the well- 
developed techniques of trapping and cooling [s, 6], and 
their degrees of freedom can also be accurately manipu- 
lated by applying laser and magnetic fields. In addition 
to new perspectives for the study of many-body quan- 
tum physics 0, 13 J these progresses offer new possibili- 
ties of atomic interferometers of unprecedented sensitiv- 
ity [9]. In particular, nonlinear dynamics, quantum en- 
tanglement and quantum metrology meet together in sys- 
tems of Bose condensed atoms [10l - [l5| , in which nonlin- 
ear dynamics may generate quantum entanglement and 
quantum entanglement may be used to implement high- 
precision quantum metrology. 

This article is a review on recent progresses in theo- 
retical studies in nonlinear quantum interferometry with 
Bose-Einstein condensed atoms. This section is an intro- 
duction for some fundamental conceptions for quantum 
interferometry, atomic BECs and nonlinear dynamics in 
BECs. In the next section, we present a description for 
MQC of Bose condensed atoms and discuss how to imple- 
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ment quantum interferometry via atomic matter- waves, 
in particular, the Bose- Josephson junctions. In the third 
section, we show how to reahze many-body quantum in- 
terferometry via an ensemble of condensed atoms. The 
last section summarizes this article and discusses some 
perspectives in this field. 



A. Quantum interferometry 

In an optical interferometer, two or more light beams 
are coherently combined for interference and then one 
can extract their relative phase from their interference 
patterns. Optical interferometers have been extensively 
used for precision measurements, surface diagnostics, as- 
trophysics, seismology, etc. There are several configu- 
rations for optical interferometers, for an example, an 
optical Mach-Zehnder interferometer splits a laser beam 
into two beams propagating along two paths and then re- 
combines the two beams for interference and extracting 
their relative phase. 

Beyond the conventional interferometry via classical 
waves, the quantum interferometry uses the wave nature 
of particles to achieve higher precision limit which can 
not be reached by classical interferometry. For an exam- 
ple, in a quantum Mach-Zehnder interferometer, a quan- 
tum state is transferred into a coherent superposition of 
two states by an unitary operation, then the two states 
accumulate different phases in a free time evolution, and 
lastly the two states are coherently recombined for in- 
terference and extracting their relative phase. The quan- 
tum interferometry has demonstrated by various systems 
from individual microscopic particles (photons, neutrons, 
electrons, atoms, molecules, etc) to macroscopic objects 
of an ensemble of microscopic particles (superfluids, su- 
perconductors, atomic BECs, etc). 

It has demonstrated that quantum states in particu- 
lar entangled states offer unprecedented advantages for 
metrology [16"25!]. The quantum states have been ex- 
tensively used in high-precision measurements, such as, 
quantum states of photons have been used for high- 
resolution imaging, quantum states of spins have been 
used for high-resolution measurements of magnetic field, 
and quantum states of atoms have been used to build 
clocks of ultimate accuracy. Moreover, multi-atom sys- 
tems in particular entangled multi-atom systems, such 
as ultracold atoms in optical lattices and quantum 
atomic gases, provide new opportunity for designing 
next-generation atomic clocks |26j . 



B. Atomic Bose-Einstein condensates and 
nonlinear dynamics 

In an atomic BEG, the matter- waves of individual atoms 
overlap each other and therefore all atoms lose their 
distinguishability. The indistinguishability of individ- 
ual atoms means the appearance of MQC. Under such 



a low temperature, the atom-atom interaction in a BEG 
is dominated by the s-wave scattering. Under the mean- 
field description [H-Q , an atomic BEG obeys a nonlinear 
Schrodinger equation, the Gross-Pitaevskii equation, in 
which the nonlinear strength is proportional to the s-wave 
scattering length. Due to the intrinsic nonlinearity from 
the atom-atom interaction, the de Broglie matter waves 
of Bose condensed atoms have led to the development of 
nonlinear and quantum optics with atoms [27| , which are 
analogues of conventional nonlinear and quantum optics 
with photons. 

It has demonstrated the existence of many nonlin- 
ear phenomena such as four-wave mixing, solitons and 
chaos in atomic BEGs. In an optical four-wave mixing, 
a fourth light wave is produced by mixing three different 
light waves under the condition of energy and momentum 
conservation. Similarly, four-wave mixing process may 
take place in atomic BEGs [l^, in which three different 
matter waves of atoms mix to produce a fourth mat- 
ter wave of atoms. Solitons are stable solutions which 
maintain their shape unchanged during their propaga- 
tion. Due to the balance between dispersion and non- 
linear atom-atom interaction, matter wave solitons may 
appear in atomic BEGs. The dark [2^ [s^l and bright 
solitons [31I, d^l have been observed in atomic BEGs with 
repulsive and attractive atom-atom interactions, respec- 
tively. Due to the intrinsic nonlinearity from the mean- 
field interaction, macroscopic quantum chaos has been 
theoretically predicted in time driven systems l33l . |34| 
and high-dimensional systems of atomic BEGs [Sa, [Sg . 
The strength of nonlinearity (i.e. the s-wave scattering 
length) can be tuned by Feshbach resonances |37| and so 
that the nonlinear dynamics can be modulated (38l [39| . 



C. Quantum interferometry via Bose condensed 
atoms 



Since the first realization of Bose-Einstein condensation 
in ultracold atomic gases, the condensed atoms have been 
extensively used for performing quantum interferome- 
try l40l - |42j . There are two fundamental schemes for 
quantum interferometry with condensed atoms: space- 
domain interferometry and time-domain interferometry. 
In the space-domain scheme, two or more atomic BEGs at 
different spatial positions are released for expansion and 
then the matter waves from different condensates overlap 
each other. In the time-domain scheme, by means of a 
double separated oscillator technique, the Ramsey oscil- 
lations of condensed atoms involving two or more internal 
hyperfine levels have been demonstrated. 
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II. MACROSCOPIC QUANTUM COHERENCE 
AND ATOMIC MATTER- WAVE 
INTERFEROMETRY 

A. Macroscopic quantum coherence and mean-field 
description 

It is well known that wave-like properties of single par- 
ticles such as photons, electrons and atoms have been 
demonstrated in various experiments. Fantastically, for 
an atomic BEC, the collection of atoms may also be- 
have as an entire object of wave-like properties [H-Il]- 
These wave-like properties among an ensemble of par- 
ticles are a signature of macroscopic quantum coherence 
(MQC), which relates to the superposition of multiple 
macroscopic quantum states. 

In the quantum field theory, an ensemble of ultracold 
Bose atoms confined within an external potential obeys 
the many-body Hamiltonian l-A i43ll4^. 



H 



2m 



+ \J dr(ir'*t(r)^t(i.')y(r _ r')^(r')*(r), (1) 

with and ^ denoting bosonic fields of atoms. Here, m 
is the atomic mass, V^xt stands for the external potential, 
and V{r—r') describes the interaction between two atoms 
at positions r and r'. 

The time evolution of the field operators is given by 
the Schrodinger equation. 



2to 



+ V,xtir) 



dr'i/\r',t)Vir~r')^{r',t) 



*(r,i). (2) 



Solving the above equation for an ensemble of particles 
involves heavy numerical work. For an atomic BEC dom- 
inated by condensed atoms, mean-field (MF) theory al- 
lows one to well understand most behaviors related to 
MQC. However, as will be discussed in Section 3, the 
theoretical treatment must go beyond MF and take into 
account quantum fluctuations if non-classical states such 
as squeezed states and macroscopic coherent superposi- 
tions are encountered. 

The field operator ^{r,t) can be expressed in form of 



*(r,0 =$(r,t) + «''(r,i). 



(3) 



where $(r, i) = (5'(r,i)) is called as the order param- 
eter or wave-function of the condensate" and ^'(r,t) 
describes the fluctuations around the condensate state. 
Therefore, the condensate density reads as no{r,t) = 

Under ultra-low temperatures, if the average inter- 
particle spacing is sufficiently large, the atom-atom in- 
teraction is dominated by the s-wave scattering. That 



is, the atom-atom interaction is effectively described by 
a "contact interaction". 



V{r' - r) = g5{r' - r). 



(4) 



with g 



47rfi o 



and a being the s-wave scattering length. 



Physically, g > and g < correspond to repulsive 
and attractive interactions, respectively. Inserting the 
delta interaction potential into the many-body quantum 
Hamiltonian (1), one can find 



H 



2m 

dr*^(r)*t(r)^(r)^(r) 



*(r) 



(5) 



For an atomic gas dominated by condensed atoms, the 
fiuctuation part ^'(r, is negligible. 

Under conditions of ultralow temperature and dilute 
density of atoms, that is, (a) ^'(r, « ^{j^,t) of very 
small fluctuations, (b) the inter-atom interaction is a con- 
tact interaction, and (c) the average inter-particle dis- 
tance is much larger than the s-wave scattering length, 
an atomic BEC can be described by a mean-field Hamil- 
tonian, 



MF 



dr$*(r) 



+ Vextir) 



+ 7 



2m 

rfr$*(r)$*(r)$(r)$(r). 



$(r) 



(6) 



The time evolution of the condensate wave-function 
<I>(r) obeys a nonlinear Schrodinger equation, the Gross- 
Pitaevskii equation (CPE), 



2m 



V,xt{r)+g\Hv,tt 



$(r,t). 



(7) 

Here, the term of g describes the nonlinear mean-field 
interaction between atoms. For a stationary state, 
<f>(r,t) = (l){r) exp{—iij,t/h), the spatial wave function 
obeys a time-independent GPE, 



/^0(r) 



2m 



+ K:rt(r) +#^(r) 



</)r), (8) 



with fi denoting the chemical potential. 



B. Atomic matter- wave interference 

The first interference experiment of two atomic BECs was 
implemented in 1997 by Ketterle's group at MIT j45| . 
In this experiment, two atomic BECs released from a 
double-well potential freely expand and then form clear 
interference patterns. Generally speaking, the double- 
well potential can be created by superposing a poten- 
tial barrier on a harmonic potential. In MIT's experi- 
ment, the harmonic trap is a magnetic potential and the 
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barrier is provided by a blue-detuned far-off-resonance 
laser beam. There are several other schemes for forming 
double-well potentials, such as, radio- frequency-induced 
adiabatic potentials created by atom chips [ia-'is] and 
superposing a weakly optical lattices on a strongly har- 
monic trap 0, HO] ■ 

There are two different methods for preparing two 
BECs in a double-well potential. The first method is 
loading cold atoms into a double-well potential and ob- 
taining two BECs via evaporative cooling. The second 
method is preparing a BEC in a harmonic trap and then 
splitting it into two BECs via increasing the potential 
barrier. 

Interference patterns will appear if atomic matter 
waves from different condensates may reach same spatial 
positions. A natural way is switching off the external 
potential. The experimental data shows the interference 
patterns differ from the ones of two point-like monochro- 
matic sources and two point-like pulsed sources. This 
indicates that the two BECs released from a double-well 
potential can not be simply treated as two point-like 
monochromatic sources or two point-like pulsed sources. 
In general, for particles from two point-like sources, their 
fringe period of the interference pattern is the de Broglie 
wavelength associated with the relative motion of parti- 
cles, A — h/p — ht/{md), where h is Planck's constant, p 
is the relative momentum, t is the time, d is the distance 
between two sources and m is the particle mass. For the 
case of two BECs released from a double-well potential, 
the fringe period formulae for two point-like sources is 
only valid for the central fringe [45| . 

For simplicity, we consider an elongated system of con- 
densed atoms. That is, the atoms are confined by a 
strong transverse confinement, ■mujp{y'^ + z'^)/2. Integrat- 
ing the transverse coordinates, the elongated system can 
be described by a one- dimensional (ID) GPE [5l|,l52|: 

d 9 
ih—^{x,t) ^ Ho<^{x,t) + \ \^{x,t)\^ ^{x,t), (9) 

where i7o = —{fi? /2m){d'^ / dx'^) + V{x,t), A characterizes 
the nonlinear MF interaction, and V{x, t) is an external 
potential. 

Usually, the condensate wave-function $(r,t) is nor- 
malized to the total number of condensed atoms N . To 
compare with the conventional Schrodinger equation of 
no MF interaction, we normalize ^{x, t) for our ID model 
to one and its nonlinearity strength A = 2NasUjph is de- 
termined by iV, the s-wave scattering length as and the 
transverse trapping frequency Wp [S^l. For convenience, 
we introduce the dimensionless model equation by choos- 
ing the natural units of m = ?i = 1. The external po- 
tential V{x^t) is formed by a superposition of a time- 
independent harmonic potential and a time-dependent 
Gaussian barrier, see Fig. 1. That is, 

y(a;,t) = ic.V + S(t).exp(^-^), (10) 

where lo is the trapping frequency, d is the barrier width. 



and the barrier height B(t) is a function of time. If 
B{t) > Lo'^d'^, V{x,t) is a double-well potential with two 
minima at a; = ±dy/2\n[B/{d'^uj'^)]. In the splitting pro- 
cess, the barrier height is adiabatically increased and 
therefore the single condensate in a harmonic trap splits 
into two condensates in a deep double-well potential. To 
observe the interference pattern, one has to switch off 
the potential barrier at least and let the two condensates 
overlap each other. 




FIG. 1: The splitting process from a single condensate to two 
condensates, (a) The initial state (green soUd curve) at time 
t = in a harmonic trap (blue dashed curve) ; (b) the density 
evolution |$(x,t)|'^ for the splitting process; and (c) the final 
state (green solid curve) at time t — 100 in a double-well 
potential (blue dashed curve). 

In Fig. 1, we show our numerical simulation for the 
splitting process. Initially, the condensate stays in the 
ground state for the harmonic trap. Then the system is 
transferred into the ground state for a deep double-well 
potential by slowly increasing the barrier height. If the 
barrier is sufficiently high, the condensates in two wells 
almost have no overlap and so that they can be looked 
as two independent condensates. In our calculation, we 
choose the interaction strength A ~ 20, the barrier width 
d = \/2, the trap frequency uj = 0.2tt, and the barrier 
height B is linearly increased from zero to 25. 

The external potential V{x,t) is suddenly switched off 
when two condensates have been prepared in the double- 
potential. Thus the two condensates freely expand and 
interference fringes gradually appear when atoms from 
different condensates reach the same position. In Fig. 2, 
we show our numerical simulation for how the interfer- 
ence fringes appear. In our calculation, we choose the 
final state of the splitting process shown in Fig. 1 as the 
initial state before switching off the external potential. 
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FIG. 2: Interference fringes of two freely expanding conden- 
sates, (a) The density evolution |<l>(x,t)|^ after switching off 
the external potential V{x,t); and (b) the interference fringe 
at time t = 3. 



C. Nonlinear excitations in atomic matter- wave 
interference 

The MF theory for atoniic BECs predicts that solitons 
may exist in ID BEC systems. If the s-wave scattering 
length is negative, the atom-atom interaction is attrac- 
tive and acts as the role of self-focusing in nonlinear op- 
tics, solitons appear as shape-maintaining wavepackets 
in propagation. If the s-wave scattering length is posi- 
tive, the atom-atom interaction is repulsive and acts as 
the role of self-defocusing in nonlinear optics, solitons 
appear as shape-maintaining notches in propagation. 

In matter-wave interference experiments of atomic 
BECs, due to the intrinsic atom-atom interaction, the 
interference patterns differ from the ones of linear waves 
and nonlinear excitations may gradually appear. Theo- 
retically, it has demonstrated that dark solitons can be 
generated by collision (S^HJl and quantum tunneling [sl] 
in double-well BEC interferometers. Additionally, non- 
linear effects in the BEC interference are studied by using 
exact solutions of the one-dimensional GPE 56] . Exper- 
imentally, for atomic BECs with repulsive MF interac- 
tion, dark solitons are gradually generated by merging 
and splitting BECs [57-60]. 

The soliton generation depends on the initial phase 
difference and how fast the recombination occurs. In a 
sufhciently slow recombination, if the initial state is the 
ground state for the double- well potential of no phase dif- 
ference, the condensate is adiabatically transferred into 
its ground state for the final harmonic trap. However, if 
the initial state is the first-excited state for the double- 
well potential of tt phase difference, the condensate will 



be adiabatically transferred into its first-excited state for 
the final harmonic trap. The notch of this excited state 
can be looked as a dark soliton. Usually, more number of 
solitons will be generated in a faster recombination . 
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Time 

FIG. 3: In-trap interference and nonlinear excitations. The 
density evolutions |$(2;,t)|^ for A = (linear system), A = 
10, A = 20 and A = 40 are shown in rows (a), (b), (c) and 
(d), respectively. All initial states are chosen as the ground 
states for BECs within a double-well potential of parameters 
Lu = 0.27r, d = 1 and B = 40. 

Based upon the ID GPE (7), we explore how dark 
solitons appear via numerical integration, see Fig. 3. 
The initial states are set as the condensate ground states 
for the external potential V{x,t = 0) with parameters 
uj = 0.27r, d = 1 and B (t = 0) = 40. Following the 
experiments on recombination from a double- well poten- 
tial into a harmonic potential [stI [ssj , we only switch 
off the potential barrier and keep the trap frequency 
w unchanged. There are several schemes for switching 
off the potential barrier. In our simulation, similar to 
the matter-wave interference in free space, we suddenly 
switch off the potential barrier, that is, B {t > 0) = 0). 
Due to the harmonic forces, two condensates collide again 
and again. The interference patterns strongly depend on 
the interaction strength. For a linear system (A — 0) 
of no inter-particle interaction, the interference patterns 
periodically oscillate with the period for the harmonic 
trap. For the nonlinear systems of repulsive interactions 
(A > 0), the periodicity for the interference patterns is 
destroyed by the nonlinear inter-particle interactions and 
dark solitons are gradually generated. In general, soli- 
tons gradually appear from the overlapped regions and 
more solitons are generated in more strongly nonlinear 
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systems. 



D. Bose-Josephson junction (BJJ) 

In 1962, Brian Josephson predicted the existence of su- 
percurrent through an insulator barrier between two su- 
perconductors 61,[63]- This tunneUing effect of electrons 
is named as the Josephson effect. The device of two su- 
perconductors linked by an insulator barrier is known as 
a superconductor Josephson junction (SJJ), which has 
been extensively used to build quantum interference de- 
vices for sensitive measurements. 

To detect and exploit the MQC among BECs of atoms, 
similar to the SJJs, it is natural to link different con- 
densates with Josephson couplings. The peculiar tun- 
nelling phenomena in Josephson coupled BECs are called 
as Bose-Josephson effects, and the physical systems of 
two BECs linked by Josephson couplings are called as 
Bose-Josephson junctions (BJJs). There are two differ- 
ent types of BJJs: the external and internal BJJs 44], see 
their schematic diagrams in Fig. 4. An external BJJ [AQh 
[50| involves two condensates in a double-well potential 
and the Josephson coupling is provided by the quan- 
tum tunnelling through the barrier between two wells, 
see panel (a) of Fig. 4. An internal BJJ [63l - [65| involves 
a two-component condensate of atoms occupying two hy- 
perfine levels which are coupled by the external fields, see 
panel (b) of Fig. 4. Besides the Josephson effects in two- 
mode systems, non-Abelian Josephson effects have been 

66[ . Further- 



predicted in a spin-2 BEC of five modes 
more, the Josephson oscillations have been demonstrated 
in an array of BECs in optical lattices jG^-ffiT 




FIG. 4: Schematic diagrams for Bose-Josephson junctions: 
(a) an external Bose-Josephson junction linked by quantum 
tunnelling, and (b) an internal Bose-Josephson junction via a 
two-component BEC linked by Raman fields. 

Because of the nonlinear MF interaction, several ex- 
otic macroscopic quantum phenomena have been demon- 
strated, such as macroscopic quantum self-trapping 
(MQST) [toI - It^ and macroscopic quantum chaos [sj, 
[341 • Moreover, in contrast to the linear Landau-Zener 
(LZ) tunneling in two-state single-particle quantum sys- 
tems, loop structures may appear in the nonlinear LZ 
tunneling in two-mode BEC systems [tsI - IttI . Actually, 
these loop structures and MQST are two different sides of 
the bistability induced by the nonlinear MF interaction. 



In this subsection, under the MF description, we will 
focus on discussing the macroscopic quantum phenomena 
in BJJs. At first, we introduce an unified model for BJJs. 
Then, we discuss the dynamics in time-independent sys- 
tems, such as, Rabi oscillation, MQST and static bifurca- 
tion. At last, we analyze the dynamics in time-dependent 
systems, such as, chaos and slow passages across the crit- 
ical point. 



1. Two-mode approximation 

In the MF theory, BJJs can be treated as multiple cou- 
pled modes of macroscopic matter waves [ssl . u^ . For 
an external BJJ with a sufficiently high potential bar- 
rier, the condensate wave-function ^(r,t) can expanded 
by two Wannier states for the two wells [71,, X9]. That is. 



$(r,t) = Vi(i)0i(r)+V2(i)02(r), 



(11) 



where the spatial functions (pj (r) (j = 1 and 2) are the two 
Wannier states and the time-dependent amplitudes ipj (t) 
are complex numbers. As the condensate wave-function 
^{r,t) is normalized to the total number of condensed 
atoms N and each Wannier state is normalized to one, 
we have \Mt)f + \Mt)\^=N. 

Under the two-mode approximation (11), integrating 
all spatial coordinates, the MF Hamiltonian read as. 



H 



MF 



- J(Vi*V2 + ^2*V'l) + £1 IV'll' + £2 |V'2|' 

'Iun\^i\* + lu22\^p2f 



+ 2Ui2\^l\^\tp2\^ 



'-C/ll<t^22 



iri^2f 



1«*12 



l^ll (^^1*^2 



+ t^22«*12 \tp2\ (V'l'/'2 

with the inter-mode tunnelling strength 



^2^1) 



(12) 



J 



dr 



0;(r)(- 



2rn 



+ V)Mr) 



the zero-point energy for atoms in the j-th mode. 



Ej ^ I dr 



2m 



the intra-mode interaction strength. 



Ujj^g / dr |(/.,(r)| 



the inter-mode interaction strength. 



C/12 = 5 / dr 101 (r) I' (r 
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and the inter-mode exchange colhsion interaction 
strength, 

Ujk^im ^9 J dr [0*(r)(/)^.(r)0,(r)0„j(r)] . 

In the exchange collision denoted by jk ^ Im, an atom 
from the j-th mode collides with an atom from the k-th 
mode and then these two atoms are transferred into an 
atom in the 1-th mode and an atom in the m-th mode. 
For a sufficiently deep double- well potential, two Wannier 
states are well localized at the two well-centers, so that 
the high-order overlaps between two Wannier states are 
very small. This indicates that the terms in last four 
rows of Hamiltonian (12) can be ignored, that is, 

Hmf - WV^2 +i^;V^i) +ei iV^il' +£2 |V'2|' 

+^C/ii|V'i|Vic/22|V'2|'. (13) 

For an internal BJJ, in which a two-component EEC 
involves atoms in two-hyperfine levels coupled by laser or 
radio-frequency fields, its Hamiltonian in quantum field 
theory reads as 



H = Hi+H2 + Hi 



(14) 



with the single component Hamiltonian for the first com- 
ponent, 

Hi= J ^{r,t){-^ + V{r,t))^,{r,t)dr 

+ ^ [ i>\{v,t)i>\{y,t)^',{v,t)i>^{v,t)dv 



A 



*}(r,t)*i(r,t)dr, 



the single component Hamiltonian for the second com- 
ponent. 



H. = 



2m 



+ 



+ j J *kr,i)*2(r,t)dr, 

and the inter-component interaction and the linear cou- 
pling between two components, 

Hi2= 912 I ^[r,t)■^\[r,t)■^l{v,t)i!2{r,t)dT 



(*}(r, t)*2(r, t) + %{v, t)*i(r, t))dv 



Ml 



transition between the two hyperfine levels and J7 is the 
Rabi frequency. The symbols >i'j(r, t) and ^j{r,t) are 
Bose creation and annihilation operators for atoms in 
respectively. 

The system can be described the MF theory if the con- 
densed atoms dominate the whole system, i.e., ^j(r, t) « 
$j(r, t). Under a sufficiently strong confinement, only 
the ground state for the external potential will be occu- 
pied and so that the condensate wave-functions can be 
expressed as [13, Isif 



$,(r,t) =7A,(i)0(r) 



(15) 



with (/)(r) describes the spatial distribution and the 
complex amplitudes satisfy the normalization condition, 
\il^i{t)\^ + \tp2(t)\^ = N. Here, TV is the total number 
of atoms. Inserting the approximation ansatz (15) into 
the Hamiltonian (14) and integrating all spatial coordi- 
nates, one will find that the complex amplitudes obey the 
following MF Hamiltonian, 



Hmf = -J (V'i>2 + V^a^i) + ei IV'il + £2 IV'2 

+ ^f/ll|^l|Vic/22|V'2|^ 
+ C/l2|V'l|'l^2|', 



(16) 



with the Josephson coupling strength, 

the single-particle energy for the first-component. 



ei 



dr 



0*(r)(^ 



2m 



V 



y)'^(r) 



the single-particle energy for the second-component. 



£2 



dr 



ft2y2 ^ 

2m 2 



the intra-component interaction strength. 



9rj / dv |0(r)| 



and the inter-component interaction strength. 



U12 = 512 / dr 



I0(r)r 



The time evolution of BJJs obeys two coupled differ- 
ential equations 



- dip J (t) _ dHMF 
' dt ~ di:*Aty 



(17) 



Here, m is the single-atom mass, gjk = Airli^ajk/m with 
Gjk denoting the s-wave scattering length between atoms 
in hyperfine states \j) and \k), V{Y,t) is the external 
harmonic potential, A is the detuning to the resonant 



with Hmf given by the Hamiltonian (13) or (16). Obvi- 
ously, the time evolution conserves the total number of 
atoms, N = 1-01(^)1^ -t- 1-02(^)1^. This means that TV com- 
mutates with the Hamiltonian. Thus, one can omit the 
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terms 0{N) and 0{N'^) in Hamiltonians (13) and (16) 
and rewrite them into an unified form, 

6 E 
H = - {n2-ni) + {n2 - riif - J (-0i'(/'2 + '02'0i) , 

(18) 

with n, = ^P*^, = l^jf , d^e2-ei+N {U22 - Un) /4, 
and Ec = Un + U22 for external BJJs or _Ec = C^ii + 
U22 — 2J7i2 for internal ones. Correspondingly, the time 
evolution of BJJs obeys the following two coupled dis- 
crete GPEs, 

= -^^1 + ^ (l^iP - \U) V'l - ^^2, (19) 



ih 



dip2 
dt 



:'02 



Er 



\i'i\^)i^2-Ji'i- (20) 



Under the transformation, Ec — > —E^ J — > —J and 
5 — > —5, the Haniiltonian (18) H becomes as —H. This 
means that the diniensionless Hamiltonian Hd — H/J is 
invariant under the transformation. In this article, the 
coupling strength J is assumed to be non-negative. 

Using the atomic numbers nj and the phases 9j, the 
complex amplitudes can be expressed in the form of 
V'iCO — y^^^Pi^^j)- Thus, one can introduce two new 
variables, the fractional population imbalance 



ni - n2 ni - 712 



N 



(21) 



ni + n2 
and the relative phase 

02-91. (22) 
From Eqs. (19) and (20), one can find {0,z) obeying 
d6 ErNz 2Jz 



dt 



: COS — S, 



dz 
'di 



-2JVi - z^i 



(23) 



(24) 



where we have assumed the Planck constant h ~ 1. 
Clearly, the equations of motion for {9, z) are classical 
Hamiltonian equations for 



Hpi9,z) 



ErN 



z^ + 5z~2Jy/l- z^ 



(25) 



which describes a non-rigid pendulum. Here, (9, z) are 
a pair of canonical coordinates corresponding to the an- 
gular momentum and the angular variable for a classical 
pendulum. Thus, the unified MF Hamiltonian for BJJs 
is equivalent to a classical non-rigid pendulum Hamilto- 
nian. 

Although Eqs. (19) and (20) are equivalent to Eqs. 
(23) and (24), numerical divergence appears in Eq. (23) 
if z = ±1 [i.e. all atoms stay in only one of the two 
modes]. Therefore, for systems involving 2 = ±1, one has 
to choose Eqs. (19) and (20) for numerical simulation. 



2. Rabi oscillation, macroscopic quantum self-trapping and 
plasma oscillation 

The stationary states of BJJs described by Hamiltonian 
(18) correspond to the fixed points for the classical non- 
rigid pendulum described by Hamiltonian (25). Math- 
ematically, the fixed points for (9, z) are determined by 
^ = and ^ = 0, and their general solutions can be 
written in forms of Jacobian and Weierstrassian elliptic 
functions (zUll- 

Similar to the electronic current in a SJJ, one can de- 
fine the atomic current in a BJJ as, 



d{n2 - rii) 
Jt 



dz 



= N— = -2NJ\/l - z2 sin 61. (26) 



For the BJJ described by the Hamiltonian (18), due to 
the conservation of its total number of atoms in two 
BECs, there is no phenomenon of a direct current be- 
tween two BECs. Although its atomic currents are al- 
ternating currents, its phase evolution is not as same as 
a SJJ with AC Josephson effects, in which the relative 
phase varies linearly with time. Firstly, in addition to 
the running-phase modes induced the asymmetry 6 cor- 
responding to the bias voltage in a SJJ, as mentioned 
above, running-phase modes can also appear in a sym- 
metric BJJ with strong nonlinearity. Secondly, a BJJ 
may support Josephson effects of AC atomic currents and 
oscillating relative phase. By loading atomic BECs into 
an asymmetric double-well potential and a washboard 
potential, the AC and DC Josephson effects in atomic 
BECs have been experimentally demonstrated 82] . 

For a given system of parameters TV, Ec, J and S, from 
Eqs. (23) and (24), its stationary states {9*,z*) can be 
obtained by solving 



ErNz 



2Jz 



: COS 9 — 5 = 0, 



2Jy/l 



z^ sin 61 = 0. 



(27) 



(28) 



Obviously, from Eq. (28), all stationary states request 
the relative phase 9 = kw with integers k. The states of 
even integers k = 2j and odd integers fc = (2j -I- 1) are 
called as in-phase and 7r-states, respectively. For simplic- 
ity, we will focus on discussing the region of — tt < 9 < 
+n. 

(1) Rabi oscillation. - For a symmetric system without 
nonlinearity, (5 = and Ec = 0, there are two stationary 
states {9* = 0,z* = 0) and {9* = tt, z* = 0), which have 
no population imbalance. Its general solutions are pe- 
riodic oscillations around one of two stationary states 
with a frequency loh = 2 J. Actually, this linear system 
is similar to the Rabi model for a two-level system cou- 
pled by classical lasers [11]. Therefore, these periodic 
oscillations around a stationary state are called as Rabi 
oscillations. In Fig. 5, we show the Rabi oscillations in a 
linear Bose- Josephson system described the Hamiltonian 




(18) oi J ^ 1, S ^ 0, Ec ^ and N ^ 1000. In which, 
the black solid dots denote the stationary states. 

(2) MQST. - For a symmetric system with non-zero 
Ec, in addition to the two normal stationary states with 
no population imbalance, two new stationary states ap- 
pear if E^N^ > 16J^. If Ec > 0, two new stationary 

states are 7r-phase states, (j)* — Tr,z* — ±y^l — -gl^j . 

If Ec < 0, two new stationary states are in-phase states, 

(^9* = 0, z* = ±y^l - -gl^j • The non-zero population 

imbalance is a direct signature of MQST (toI - It^ . The 



critical point, where the two new stationary states appear 
in addition to the normal stationary state, is as same as 
the bifurcation point for a Hopf bifurcation 'si'] . Usually, 
the general solutions are periodic oscillations enclosing at 
least a stationary state. However, due to the nonlinear- 
ity, there may exist a new type of solution, running-phase 
mode, in which the population oscillates periodically and 
the relative phase changes monotonously. 

In Fig. 6, we show the Rabi oscillations and MQST 
in a nonlinear Bose-Josephson system described by the 
Hamihonian (18) oi J ^ 1, 6 = 0, Ec = 0.01 and 
N — 1000. The stable and unstable stationary states 
are denoted by black and red solid dots, respectively. In 
addition to the normal stationary states of no popula- 
tion imbalance, there appear nontrivial stationary states 
of non-zero population imbalance, which are self-trapped 
states. Correspondingly, there are three different types of 
non-stationary states: (i) Rabi oscillations in the green 
orbit of {z)t = and {d9/dt)t = 0; (ii) running-phase 
MQST between the green and red orbits of {z)t ^ and 
{de/dt)t ^ 0, which is labeled as MQST-A ; and (ui) 
MQST in the red orbits of (x)t ^ and {de/dt)t = 0, 
which is labeled as MQST-B. Here, {A)t denotes the time 
averaged value for the variable A. 



(3) Plasma oscillation. - Around a stationary state, a 
BJJ may have small harmonic oscillations corresponding 
to plasma oscillations in a SJJ 71, 72, 84]. The solution 
around a stationary state {9* , z*) can be written in form 
of 



z = z + z. 



(29) 



where z and 9 are small oscillations around {9*,z*). In- 
serting this solution into Eqs. (23) and (24) and keeping 
only the linear terms of z and 9, one will obtain the fol- 
lowing linearized equations. 



d9 _ EcN „ 2 J cos ( 



3/2 ' 



2Jsin6'* 

l-iz*f 



dz 2Jz*sm9* 



dt 



2Ja/1 



■ cos 9* 



I (30) 



(31) 



As all stationary states satisfy sin 6** — 0, the terms of 
sin0* disappear in the linearized equations. 

For a linear system, due to i?c = and z* = 0, the 
linearized equations read as 



d9 
'dt 



= [2Jcos9*]z, 



dz 
'dt 



= -[2Jcos9*]9. (32) 



Clearly, these small oscillations are Rabi-like oscillations 
of an angular frequency ujp — 2 J, which is just the Rabi 
frequency. 

For a nonlinear system around the normal stationary 
states of z* = 0, the linearized equations become as 



d9 
'di 



EcN 



2 J cos 61* 



z, [2Jcos9*]9, (33) 
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which describe the small oscillations of an angular fre- 
quency ojp — 2 1 + 4j^o^g. ■ Mathematically, for sta- 
tionary states satisfying ^f^^g, < — 1, the angular fre- 
quency uip may take imaginary values. The appearance 
of imaginary ujp means instability of a stationary state 
and divergence of all small perturbations to this station- 
ary state. Thus, for an unstable stationary state, there 
is no small oscillation around it. 

For a nonlinear system of -Ec > and (5 = 0, the in- 
phase stationary state of cos9* — +1 and z* = is sta- 
ble, but the anti-phase stationary state of cos^* = — 1 
and z* = becomes unstable if |-Ec^| > 4J. For a 
nonlinear system of i?c < and S — 0, the anti-phase 
stationary state of cos9* — —1 and z* = is stable, but 
the in-phase stationary state of cos 6** = +1 and z* = 
becomes unstable if |-Ec-^| > 4 J. Actually, the transition 
from stable to unstable of a normal stationary state ac- 
companies with the appearance of two additional stable 
stationary states of non-zero population imbalance. This 
confirms that the bifurcation from normal to self-trapped 
states is a Hopf bifurcation 81]. 

For a strongly nonlinear system of |i?c-^| > 4 J and S = 
0, the small oscillations around the stationary states of 

non-zero population imbalance z* = ±wl 



WW obey 



de 



Er-N 



16J2 



1 



dz 
'dt 



8J 



2 1 



ErN 



(34) 



These small oscillations are sinusoidal oscillations of an 
angular frequency ujp — 2J\ 



16J2 



- 1. 



Different from a linear system, the plasma frequency 
of a nonlinear system not only depends on the Joseph- 
son coupling J, but also depends on the nonlinearity Ec 
and the relative phase 9*. In a strongly nonlinear sys- 
tem, there are phase-dependent stationary states of self- 
trapped non-zero population imbalance and small oscil- 
lations around these stationary states. Due to the com- 
petition between the nonlinearity and the Josephson cou- 
pling, the plasma frequency gradually decreases to zero 
when the system approaches to the critical point, where 
self-trapping appears. 



3. Shapiro resonance and chaos 

In the presence of a periodic driving field, if the energy 
difference of a particle (or a composite particle such as a 
Cooper pair) at two sides of the potential barrier can be 
compensated by the driving field, the system exhibits res- 
onant tunnelling phenomena named as the Shapiro res- 
onance [t^, [13] ■ In a SJJ subjected to a bias voltage, 
V = Vo-\-Vi cos [uJdt) , Shapiro resonances occur when the 
driving frequency satisfies nujd — uJac- Here, n are posi- 
tive integers and uJac = ^eVg/h is the angular frequency 
for the AC Josephson oscillations in a SJJ with a DC 
bias voltage V = Vq. The resonance condition indicates 
that the energy difference is compensated by the energy 



of n "photons" from the driving field. This means that 
the physical pictures for Shapiro resonances and photon- 
assisted tunnelling [85l - l87| have no any difference. 

In a BJJ, the asymmetry 6 acts as the role of the bias 
voltage in a SJJ. Similarly, we assume the asymme- 
try has a constant component and a periodic component, 
that is, S = Sq + Si cos LUdt. Such an asymmetry can be re- 
alized by controlling double- well asymmetry (or detuning 
to resonance) for external (or internal) BJJs. Introducing 
a transformation 



(fii (t) = ipi (t) ■ exp 



¥^2 (t) = 1p2 {t) ■ cxp 



cos {uJdt) dt 



iSi f , , , 
— Y \ cos [iOdt)dt 



(35) 



(36) 



the coupled discrete GPEs (19) and (20) read as 



= +1^2 + ^ {W2? - l^iP) ^2 - JVi, (38) 
with the new coupling strengthes 



J = J exp 
= J cxp 



+i5i I cos {iOdt) dt 



-i— sin (ujdt) 



and J* is the complex conjugate for J. Using the ordi- 
nary Bessel functions (x), the coupling strength J can 
be expanded as. 



(39) 



with X — Si/ujd- Thus the Hamiltonian (18) becomes 
H = ^ (n2 - ni) + ^ (^2 - nif 
- J J2 Bk{x){e'''^-'iplip2+e-'''^-'ip;^i),m 

fc— — oo 

with rij = fjfj = ■ 

For a linear system of zero Ec, if coupling term can 
be regarded as a small perturbation, Shapiro resonances 
are expected to take place when the energy of an inte- 
ger number of "photons" matches the static asymmetry, 
nhujd = So- Near a n-photon resonance, the detuning 
A (uj) — nujd — S is very small. Using the rotating-wave 
approximation, in the vicinity of a n-photon resonance, 
the fast oscillating terms with fc 7^ n in the above Hamil- 
tonian can be ignored. That is, the inter-mode coupling 
is dominated by the term of _B„ {5i /uJd)- Therefore, at the 
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n-photon resonance, the driven asymmetric system is ap- 
proximately equivalent to a undriven symmetric system 
of an effective coupling strength J^g = JBn {5i/ujd). 

For an actual system under a weak driving, a resonance 
may occur if a rational multiple of the driving frequency 
matches the intrinsic frequency of the undriven system. 
Since the intrinsic frequency depends on both the cou- 
pling strength and the nonlinearity strength, the actual 
resonance positions are shifted away from nfiiod = [13 ■ 
Generally, if resonance condition is not satisfied and the 
driving is sufficiently weak, the oscillations are quasi- 
periodic ones including two frequency component: the 
intrinsic oscillation frequency and the driving frequency. 

We now consider systems of strong driving fields which 
can not treated as perturbations. Under a strong driv- 
ing, the time evolution of a linear system will be quasi- 
periodic oscillations, which include the intrinsic oscilla- 
tion frequency and the driving frequency [88| . For a non- 
linear system, chaotic oscillations may appear if the driv- 
ing is sufficiently strong [H, [sl, HI] . In Fig. 7, we show 
Poincare sections for a driven nonlinear Bose-Josephson 
system (18) with J ^ I, = 0.01, N = 1000 and 
5{t) = 5i cos (27rt). In our simulation, for a given initial 
condition, we record a data per driving period. In the 
Poincare sections, a curve or a circle of infinite dots cor- 
responds to a regular motion and a region of infinite dots 
corresponds to chaotic motions. For a weak driving, the 
oscillations are regular, see panel (a). When the driving 
strength 5i increases, chaos appear in the vicinity of the 
unstable stationary state for the undriven system and 
the region of chaotic motions increases with the driving 
strength. However, even in the sea of chaos, it is possible 
to find some islands of regular motions. In our figure, 
the chaotic sea is denoted by the region of green dots, 
regular motions in the small regular islands surrounded 
by chaotic sea are denoted by curves/circles of red dots, 
and other regular motions are denoted by curves/circles 
of blue dots. 

It is also possible to find exact Floquet states for sys- 
tems driven by some periodic modulations j90l - [9^ . These 
Floquet states have no linear counterparts due to their 
existence requires a nonzero MF interaction ^Qj. These 
nonlinear Floquet states may be stable or unstable. Usu- 
ally, the stale Floquet states stay around the node points 
of the undriven system and the unstable Floquet states 
are enclosed by the chaotic motions around the saddle 
points of the undriven system. There is a triangular 
structure in the quasienergy band, which corresponds to 
the onset of a localization phenomenon called as the co- 
herent destruction of tunneling .91.] . 



4. Universal mechanism for slow passages through a critical 

point 

Mathematically, slow passages through a bifurcation 
point exhibit memory effects induced by bifurcation de- 
lay [9^ [13] • Physically, slow passages through a critical 




(b) 5., =1/8 




(c) 5, = 1/4 



(d) 6. = 1/2 





FIG. 7: Poincare sections for a Bose-Josephson system (18) 
with J = 1, Ec = 0.01, A'^ = 1000 and time-dependent asym- 
metry S{t) = 5iCos(27rt). (a) = 0.001, (b) Si = 1/8, (c) 
(5i = 1/4 and (d) Si = 0.5. 



point of vanished excitation gap connect with the break- 
down of adiabaticity and the appearance of non-adiabatic 
excitations (also called as defect modes) [IH]- When the 
system approaches to the critical point, the adiabatic- 
ity breaks down and the generated defect modes follows 
Kibble-Zurek (KZ) mechanism [osj - losj . As mentioned 
in Section 2.4.2, for a BJJ described by the Hamilto- 
nian (18), there is a phase transition from normal to 
self-trapped states, which corresponds to a Hopf bifur- 
cation fsij. Below, we present the slow passage through 
the critical point between the normal and self-trapped 

states [ggj . 

For simplicity, we assume the Josephson coupling 
strength J > and only discuss the MF ground-states 
for symmetric BJJs in this section. For a symmetric 
BJJ of positive Ec, the normal-to-MQST transition oc- 
curs in anti-phase states, which are not its MF ground- 
states. For a symmetric BJJ of negative Ec, the normal- 
to-MQST transition occurs in in-phase states, which are 
its MF ground-states. To push the system across the 
critical point, one may adjust th e nonlinearity Ec via 
Feshbach resonances [33, llOOl . Il0l| or vary the Josephson 
coupling strength J. For external BJJs, the Josephson 
coupling strength J can be controlled via tuning the po- 
tential barrier 46 50]. For internal BJJs, the Josephson 
coupling strength J can be controlled via tuning the in- 
tensity of the coupling lasers (gsI, l63] . 

To explore the normal-to-MQST transition in ground- 
states, one has to select a system of Ec < 0, 
whose critical point is determined by EcN = — 4J. 
For this system across the critical point, its ground 
state changes from {'4'2^,''l'i^) = {\/L, \^) exp{i9j) to 
(Vi - Is, VL + Is) exp{i0j) with L = N/2 and 1^ = 
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±yfL^ -U^/El. Clearly, the normal-to-MQST transi- 
tion is a type of spontaneous s ymm etr y br eaking related 
to internal degrees of freedom [8ll . lool Il02| [ . 

In a slow passage through the critical point, there are 
two characteristic times: the reaction time and the tran- 
sition time. The reaction time, Tr{t) = h/Ag{t), de- 
scribes how fast the system follows eigenstates of its in- 
stantaneous Hamiltonian. The transition time, Tt{t) = 
Ag{t)/\dAg{t)/dt\, depicts how fast the system is driven. 
The adiabatic condition keeps valid if Tr{t) < Tt{t). Oth- 
erwise, if Tr{t) > Tt{t), the adiabatic condition becomes 
invalid and defect modes appear. 

The Bogoliubov theory is a useful method for analyz- 
in g the collective excitations over the ground states [l|- 
R Il03j |. In a BJJ described by Hamiltonian (18), due to 
its spatial degrees of freedom have been frozen, the exci- 
tations will not change the spatial distribution but only 
change the population distribution. Therefore, under the 
MF description, the perturbed ground state read as 



= [X] + S^jit)] exp ~-^lt 



(41) 



with tjjf^ = Xj exp (— and the ground state chem- 
ical potential ii. Inserting the above perturbed ground 
state into coupled GPEs (19) and (20), one will obtain 
the linearized equations for the perturbations, 

with Dirac kets 



The two coefficient matrices 



have elements 



911 



qii qi2 

921 922 



El 

4 



[P] = 



(2|xi| 



Pll Pl2 
P21 P22 



\X2\ 



(44) 



912 = - J - -^X2Xi, 



E, 



921 = - J- —X*iX2, 



922 



El 

4 



(21 



X2\ 



IXil 



E EE 

Pll = -^Xl, Pl2 = P21 = ^X2Xl, P22 = -^Xl- 

As the complex amplitudes Xj fo^ the ground state are 
determined by the parameters, given the formula above, 
all matrix elements are determined by the parameters. 



The perturbations Sijjj (t) can be written in form of 

Sipj{t) = ^ [uji exp{~iuit) + v*i exp {+iuit)] (45) 
I 

with the excitation frequencies uji, the complex ampli- 
tudes {uji,Vji) and the complex conjugate v*i for Vji. 
Substituting the expansion (45) into Eq. (42), one can 
find 



hhJi \u)i = [q] 



[P] 



(46) 



by comparing the coefficients for terms of exp(— ia;;^). 
Similarly, comparing the coefficients for terms of 
exp(-l-jwit), one will get 



-hu:i\v*),^[q]\v*), + [p\\u*),. 
Here, the Dirac kets are defined as 



Pll 
hi 



~ [ P21 



(47) 



(48) 



for P = u and v. Combine the above two equations to- 
gether, one will obtain the eigen-equation for the excita- 
tion modes. 



I - 



[9] 

-[p*\ 



[p] 

-[9*] 



(49) 



In which, [A*] is the conjugate matrix of [A\. 

Solving the eigen-equation for a symmetric system who 
has 5 = 0, one will find that, apart from a trivial gapless 
excitation mode, there exists a gapped excitation mode 
over the ground state. The energy gap for the gapped 
excitation reads as 



Ag = hWg = 



2Ja/1 



E,N 
4.7 ' 



for J > Jc, 



2J 



^1^-1, forJ<J„ 



(50) 



where Jc = \EcN /4\ — —EcN/A is the critical Josephson 
coupling strength. Clearly, the excitation gap Ag grad- 
ually decreases to zero when J approaches to Jc. Addi- 
tionally, it is easy to find that the excitation frequency 
LUg for the gapped mode is just the Plasma oscillation 
frequency cUp discussed in Section 2.4.2. This is because 
that both ujg and ujp are obtained from the linearized 
equations for the perturbations. 

We now discuss the dynamical mechanism near the 
critical point J = Jc. We assume J{t) = Jc(l ± t/rq) = 
Jc ± at is a linear function of time t, where Tq is the 
quenching time. The relative coupling, e = \[J{t) — 
Jc]/Jc\ = \t\/T„, corresponds to the relative temperature 
in KZ theory pj5l - [98| . The vanished gap at the critical 
point means the divergence of the reaction time. In a slow 
passage through the critical point, the adiabatic evolu- 
tion becomes invalid if Tr{t) > Tt{t) [95l - [98| . The slow 
passage divides into three regions: one non-adiabatic re- 
gion sandwiched by two adiabatic ones, see Fig. 8. The 
two boundaries are determined by Tr(t) = Tt{t). 
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In the normal region, J{t) > Jc, introducing j{t) 
J{t) — Jc, we have the reaction time, 



1 



1 



2^J(J + 
and the transition time, 

2J{J+Jc) ^ 

Jc[2J+Jc) 



To, 



(51) 



n 



2e{l + e)^ 
l + 2e ' 



(52) 



where tq = 1/ (2Jc). By solving = r^, one can find the 
relation, 



2e + 1 



2[£(£ + l)]3/: 



(53) 



For a slow passage, ^ 1, using the Taylor formulation, 
one will find 



>-2/3 2/3 2/3 



2-^l\y\\l\ (54) 



Similarly, in the MQST region, J(t) < Jc, introducing 
J{t) = Jc — Jit), we have the reaction time. 



Tr = 



2^J{2Jc-J) V^(2-e) 
and the transition time, 

_ J{2Jc -J) _ e{2 - e) 

- Jc{J. - J) 

Solving Tr = Tt, we have the relation 

l-£ 



To, 



(55) 



For a slow passage, ^ 1, it is easy to find 



e 



2-\ll\fl\ \t\^2 



-12/3 1/3 
To Tq ■ 



(56) 



(57) 



(58) 



Based upon the above analysis on slow passages 
through the critical point, the universal scalings e ~ 



2/3 -2/3 
^0 Tq 



and \t\ 



anism [6,8]: e 

1/(1 + 21^) ZV/(1 + ZV) 



~ Tq^^tI^'^ recover the KZ mech- 
- ^V(i+-)^-V(i+-) and \t\ ^ 

with z = 1 and v = 1/2, for t he con- 
tinuous quantum phase transitions [H, Il04l4l06{ . This 
means that the normal-to-MQST transition is a continu- 
ous quantum phase transition. 

In Fig. 8, we show our numerical results for a Bose- 
Josephson system (18) with Ec = -0.01, N = 1000 and 
(5 = 0. The initial state is chosen as a MQST ground 
state for J = 1.5 and then J is linearly swept through 
the critical point of J = Jc — 2.5. The reaction time 
only relies on the parameters and diverges at the criti- 
cal point. But the transition time Tt depends on both 
the parameters and the quenching time Tq. For a larger 
Tq, a smaller non-adiabatic region will be found. Due 



to the divergence of induced by the gapless excitation 
at the critical point, the adiabaticity is always broken 
down when the system approaches to the critical point. 
The defect excitations appear after the adiabatic condi- 
tion becomes invalid. The oscillation amplitudes decrease 
with the increase of Tq. That is, the fast passage pro- 
cesses of shorter quenching times have large oscillation 
amplitudes. 




FIG. 8: Dynamics of slow passages through the critical point 
for the normal-to-MQST transition in a Bose-Josephson sys- 
tem (18) with = -0.01, iV = 1000 and 5 = 0. The 
above panel shows how the reaction time Tt and the tran- 
sition time Tt depend on the Josephson coupling strength J 
and the quenching time Tq. The bottom panel shows the non- 
adiabatic dynamics of the fractional population imbalance z. 



In addition to the above slow passage through the 
critical point induced by varying the ratio between the 
Josephson coupling and the nonlinear interaction, there 
are several interesting works on adiabatic evolutions 
and non-adiabatic excitations in nonlinear LZ tunnel- 
ing processes induced by varying the asymmetry [tsI - 
llOSl j . For a system of strong nonlinearity, due to 
the appearance of the loop s truct ure [tsI [tHI , the adia- 
batic evolution breaks down |l07l | and the semic lassical 
and adiabatic limits become incommutable (l08j |. Qual- 
itatively, the adiabaticity of such a nonlinear evolution 
requires to vary slowly the paramete rs wi th respect to 
the Bogoliubov excitation frequencies |l07l | . 

Besides the studies on adiabatic dynamics in BJJs, 
the adiabatic dynamics in atom-molecular BECs has also 
been studied [l09l Ill0l | . In the absence of inter-particle 
collisions, the adiabatic condition has been derived with 
a linearization procedure, in which the col lectiv e modes 
obey the usual orthonormalization relation |l09l | . In con- 
trast to the linear quantum systems, if the nonlinearity 
from two-body collisions can not be ignored, one has t o 
consider the Bogoliubov collective excitation modes [llOf . 
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III. MANY-BODY QUANTUM 
INTERFEROMETRY 



Beyond the MF theory, the many-body quantum effects 
become important in systems of low dimensionality, sys- 
tems possessing a strict symmetry constraint, degener- 
acy, or strong interaction [3, . Examples include low di- 
mensional gases, condensates of spinor atoms, rapidly ro- 
tating condensates, condensates with strong atom-atom 
interaction, and ultracold atoms confined in double-well 
potentials or optical lattices. 

To explore the many-body quantum effects, one can 
employ a second quantization theory. That is, the sys- 
tems are treated as an ensemble of discrete quantized 
particles rather than continuous fields. Within the sec- 
ond quantization formalism, the systems of multiple cou- 
pled modes obey Hubbard- like Hamiltonians Illl4ll4 |. 
and several many-b ody quant um phenomena including 
quanturn squeezing [ll5l Ill6j |. quantum phase tr ansi 
tion [117], and dynamics and collapse and revival 111 



etc have been observed via manipulating Bose atoms 
in optical lattices. In particular, several many-body 
states of two-mode systems have been used to ii nple - 
ment high-precision quantum interferometry 10l - [l5i [iIqI - 
Il23l |. Moreover, it has also explored several fundamental 
noises in quantum interfero metry with one-dimensional 
Bose atomic gases [T23 - IT27l| . 

In this section, we focus on the many-body quantum 
interferometry via two-mode systems of Bose condensed 
atoms, which could be described by the unified model 
for BJJs. Firstly, we give the many-body model in sec- 
ond quantization and discuss some key characteristics of 
its ground states. Then, we show how to generate spin 
squeezing and then use it for quantum interferometry. 
Lastly, we discuss how to prepare NOON states and then 
use them for quantum interferometry. 



A. BJJ in second quantization 



In the second quantization theory, replacing the complex 
amplitudes tpj {ip*) with the bosonic operators bj {b'j'), 
the MF Hamiltonian (18) becomes a two- mode Bose- 
Hubbard model, 



H = -{n2- fii) 



E. 



- (n2 - ni) 



-J[hlhi+h\h2 



(59) 



b\bj and the total num- 



with the number operators hj - 

ber of atoms N — (N) ~ {hi + 712) ■ In the Fock basis 
{|ni,n2)}, any arbitrary many-body quantum states of 
the above Hamiltonian can be written in form of |^) = 
Sni,n2 ^"i."2 l'^i;'^2) with C„i,„2 denoting the proba- 
bility amplitudes for \ni,n2). Obviously, this Hamil- 
tonian conserves the total number of atoms N due to 

n,h] = 0. 



The many-body ground states of the quantized BJJ 
described by Hamiltonian (59) depend on its parameters 
and the total number of atoms N. If the Hamiltonian 
is dominated by the Josephson coupling term of J, the 
ground state is a SU(2) coherent state. If it is dominated 
by the nonlinear term of positive Ec, the ground state 
is a single Fock state t) ^'^^ even or a twin Fock 
state ^(1^,^) + I ^,^)) for odd N. If it is 
dominated by the nonlinear term of negative Ec, there 
are two degenerated ground states |iV, 0) and |0, TV), so 
that any superposition of these two states including the 
NOON state ^ (|iV, 0) -I- |0, N)) is also a ground state. 

A quantized BJJ of N atoms can be treated as a spin- 
Y particle in an external magnetic field. By introducing 
the spin operators, 







+ b{b2), 


Sy 




-b\b2), 




= i(^62 


- b\bi), 



the H amiltonia n 
model [T28l - [T30l |. 

H = 



(60) 



(59) is equivalent to a giant spin 



-BxSx — BzSz + DS^^ 



(61) 



with the transverse magnetic field Bx = 2 J, the longitu- 
dinal magnetic field B^ — and the axial anisotropy 
energy D — 



0.10 



0.05 




FIG. 9: Interaction blockade, resonant tunneling and energy 
spectrum for a quantized BJJ of A'^ = 10 and Ec/J = 10. 
The first row shows the coherence and the expectation value 
(Sz) versus the asymmetry strength S/Ec. The second row 
shows how En/Ec depends on S/Ec with En denoting the 
eigenenergy for the n-th eigenstate. 

In a quantized BJJ of weak Josephson coupling, the 
competition between asymmetry and nonlinear interac- 
tion wil l result i nteraction blockade and resonant tun- 
nelling The single-atom resonant tunneling 
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occurs when the asymmetry is balanced by the interac- 
tion energy, that is, 



P = 



2 a - af 
— P = ■■ — 



(68) 



with n — n2 — ni, n' = n'2 — n[ , rij = ^2 ± 1 and n'l — 
ni =F 1. Simphfying the balance condition, one can easily 
obtain 



E. 



6rT I Ec , i \ 



(63) 



withm = f = {-f ,-f + l,---,f-l}andiV = ni+n2. 
This means that there are N resonant peaks for the 
single-particle tunneling in a system of N atoms. Ac- 
tually, due to the positivity of i?c, the resonant peaks for 
the single-particle tunneling correspond to the avoided 
crossings between the ground and first-excited states. 
Between two resonant peaks, the quantum tunneling be- 
tween two modes is blockaded by the large gap origi- 
nated from the interaction energy Ec- This is a kind of 
interaction blockade similar to the Coulomb blockade of 
electrons in quantum dots. The phenomena of resonant 
tunneling and i nter action blockade have been observed 
in experiments jl34| . 



B. Spin squeezing and quantum interferometry 



Quantum squeezing |135l . ll36l | is a consequence of Heisen- 
berg's uncertainty relations. A state is squeezed when the 
fluctuations of one variable is reduced below the symmet- 
ric limit at the expense of the increased fluctuations of 
the conjugate variable. The quantum squeezed states of 
photons have provided i mportant applications in high- 
precision interferometry jl37l Il38j | and have stimulated 
extensive interests in both theoretical and experimental 
quantum optics. 

To illustrate the conception of squeezing, we consider 
a harmonic oscillator described by the Hamiltonian, 



-9 



1 



H = — 1 — muP'Sp' . 

2m 2 



(64) 



In units of H = 1, one can introduce the creation operator 



Therefore, we have the commutation relation [A", P] — 2i 
and the uncertainty relation AAAP> 1. The state of 
AA = AP = 1 is a coherent state |l39l |. Otherwise, the 
state is a squeezed one if A A < 1 or AP < 1. 

For systems represented by spin or angular momentum 
operators, such as the Hamiltonian (61), the commuta- 
tion relations for Sx.y,z read as 



[Si, Sj] — KijkSk 



(69) 



where eijk is the Levi-Civita symbol and {i,j,k) are in- 
dices for three Cartesian components along directions 
(e^;, Cy, Cz). The commutation relations leads to Heisen- 
berg's uncertainty relations 



AS, ■ ASj > - 



{Sk 



(70) 



In analogy to a harmonic oscillator, a coherent spin state 
(CSS) [21I, Il39j | is a collective state describing a collec- 
tion of spins that are perfectly aligned or polarized in 
one direction eV-, in which the fluctuations along two 



other orthogonal directions ASi — ASj 



1 

V2 



{Sk 



Otherwise, a spin state is squeezed if ASi or ASj is 



smaller than 



(Sk) 



m, UM UM- In language of 

quantum phase states, a spin--y CSS is a phase state 
for a two-mode Bose system of N particles [141]. For 
a double-well BEC system, the differences between MF 
and many-body qua ntum dynamics have been explored 
by using CSS's [l4^ . 

Generally, a spin-^ system can be regarded as a col- 
lective system of N spin-i particles. The large spin oper- 
ators can be constructed by the collective spin operators 



Sa 



2 ^ 

1=1 



(71) 



where a = (x, y, z) and are the Pauli matrices for the 
Z-th particle. Thus, the CSS can be expressed as a direct 
product of N single-particle states 



X - i\l p, 



2 V 2mw 

and the annihilation operator 



1 



(65) 



(66) 



These two operators obey a bosonic commutation rela- 
tion [a, a^] = 1. It is convenient to introduce two new 
dimensionless operators. 



A = \/2mijjx — a + a'' 



(67) 



N 



1=1 



cos - 1;), 



e''!' sin ■ 



It), 



(72) 



where |4,); and |t); are spin-down and spin- up states of 
(jf respectively. The means that all N particles point to 
the same direction {9, 0), which is called as the mean spin 
direction(MSD). The variance of the components perpen- 
dicular to MSD is simply the sum of the variances of N 
individual particles, and it is since there are no quan- 
tum correlations among particles. Spin squeezed states 
(SSS) are associated with the quantum correlations or en- 
tanglement among particles. It is possible to reduce the 
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fluctuations in one direction at the expense of increasing 
fluctuations in the other direction. 

There are several parameters for characterizing the 
spin squeezing [2l|. Below, we give the definitions for 
three squeezing parameters, which are mostly used in re- 
lated fields. 

• The first squeezing parameter is defined acco rding 
to the Heisenberg's uncertainty relation (l43j |. It is 



given as 



?ff = jOi ^ 7. 

\{S,)\ 



(73) 



With this definition, the state is squeezed if < 1 . 

• The second squeezing parameter is determined 
the r ninim um fiuctuation along a particular direc- 
tion [l4(| . It is defined as 



lin (^ASfij_^ 4 • min (^ASn 



j72 



N 



(74) 



where the spin length j = N/2, the unitary vector 
n± denotes a direction perpendicular to the MSD 
and the minimization is over all possible directions 
n±. This parameter is set to find out the most 
squeezed direction perpendicular to the MSD. If 
^1 < 1, the state is squeezed. 

• The third squeezing parameter is the ratio of phase 
fiuct uations fo r the considered state and a reference 
CSS [l43,[l4il. It is expressed as 



^2 - 



A<j) 



(A0) 



CSS 



(75) 



Similarly, the state is squeezed if < 1. 



Now, we discuss how to enhance measurement preci- 
sion by using spin squeezing. For an example, we consider 
a Ramsey interferometry of two 7r/2-pulses sandwiching a 
free evolution. The Ramsey interferometry is equivalent 
to a Mach-Zehnder (MZ) interferometry, in which the two 
7r/2-pulses act as the two beam splitters. Assuming an 
unknown phase (f> about to the z-axis is accumulated in 
the free evolution, the state of the Ramsey interferometer 
evolves according to 

\^{t))=U\^{0)), (76) 

with the propagation operator 

tl = exp exp (^-i4>Sz^ exp {-^^^xj 

= exp (^—i(j)Sy^ exp (^—inSx^ ■ (77) 

To estimate the phase (f), one may measure a observ- 
able Q, which has 0-dependent expectation values and 



variances. According to the error propagation formula, 
the variance of (f) is given by 



Acl) 



AQ 



d(Q)/dcj, 



(78) 



In the Ramsey interferometry, it is usually to measure Sz 
of the final state. The expectation value of Sz is 



S,) = 
f 



{^{t)\Sz\^{t)) 

(*(0)|;7+^,?7|*(0)) 



sin (/) ■ {Sx)t=o - cos • {Sz)t=o (79) 



and the corresponding variance is 
2 / ^ \ 2 

— sin 



2 



(aS-^) ^ = sin^ (f> ■ (^AS^'^ ^ ^ + cos2 • (aS-^) ^ 

{2<l>) -Covf Sx,Sz) (80) 
V / 1=0 



with 



Gov (^,, Sz) = i (Sjz + SzSx) . (81) 

Thus, we can estimate the unknown phase with the un- 
certainty 



A0 : 



(A5.)^ 



d{Sz)^/dq^ 



(82) 



If the input state is an CSS, the phase sensitivity is 

1 



(^'^)CSS = 



(83) 



which is the standard quantum limit (SQL) [T^. If the 
input state is a SSS, the phase sensitivity reads as 



(A^)sss - ^■ 



(84) 



This means that the phase sensitivity could be improved 
beyond the SQL by utilizing spin squeezing. 

Usually, one can generate spin squeezing by nonlinear 
inter-particle interactions. In 1993, Kitagawa and Ueda 
demonstrated that the squeezed states of S^n oc 1/N~^/^ 
(or can be generate d by the one-axis (or two- 

axis) twisting Hamiltonians jl40|. Recently, it has also 
demonstrated that planar spin sque ezed states can be 
prepared via attractive BECs [l46j |. Therefore, using 
these squeezed states, the phase sensitivity can beat the 
SQL provided by the unentangled states. For a two- 
component atomic BEC, it has proposed that the non- 
linear twisting can be accomplished by the s-wave scat- 
tering between atoms and the spin squeezing is a use- 
ful criterion for many-particle entanglement [lo!]. Re- 
cently, two experimental groups have prepared SSS's of 
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Bose condensed atoms by using one-axis twisting and 
demonstrated high-precision quantum interferometry be- 
yond the SQL by utihzing these squeezed stat es [l4 [l5l | . 
In a practical BEC system, particle loss |l47| and non- 
condensed atoms have significant effects on the optimal 
spin squeezing achievable from the inter-atom interac- 
tions pSfisaj . 

Additionally, other con trolled in teractions such as 
laser-induced interactions , co ntinuous quan- 

tum non-demolition measurements |157l | and Coulomb 



158l - [T60j are also able to generate non- 



interactions 

classical entangled states for implementing high-precision 
interferometry. 



C. NOON states and quantum interferometry 



A NOON state |161l4l63j is a many-body quantum state. 



1 



INOON) = — (|7V), |0), 



|0>jA^)j 



(85) 



which represents a superposition of N particles in mode a 
with zero particles in mode b, and vice versa. If the two 
modes are regarded as two possible paths for the parti- 
cles, the NOON state is also called as the path-entangled 
state. Due to only two modes are involved, the NOON 
state is equivalent to the N-particle GHZ state, 



I GHZ) 



V2 ( 



N 



N 



N 









(86) 



which is a maximally entangled state. 

Because of their entanglement characteristics, NOON 
and GHZ states can be used to implement Heisenberg- 
limited interferometry [l6l|. For an example, we consider 
a N-particle Ramsey interferometer. If we input a GHZ 
state, 



|GHZ), 



N 



N 




Y' 





(87) 



after the Ramsey process, the final output state is 
|GHZ)^= 




For this final state, because the initial state has zero ex- 
pectations for Sx and Sz, it is not possible to extract 
any information about </> by measuring Sz- To esti- 
mate the phase 0, one may measure the parity operator 
P — Y\f^i erf, where of is the z-component Pauli matrix 
for the 1-th particle jl64| . With some analytical calcula- 
tions, one can easily find that the final expectation values 
and variances of P after the Ramsey procedure are given 
as 



cos 



{N(j}) , ( AP) ^ ^ sin^ {Ncj}) . (89) 



Therefore, the phase sensitivity is determined by 



d{P) ^/d(t> 



TV' 



(90) 



which reaches the Heisenberg limit 1/A^. 

To estimate from the final state (88), one can also 
measure the operator A — |'f|-)(-IJ.| -I- |-IJ-)('fl"l with |f|-) = 
It'+t)^ ^'^'^ \^) ~ It'^t)j;- "^^^"^ operator has simi- 
lar expectation values and variances for the parity oper- 
ator P. 

Although, theoretically, NOON states may enhance 
the measurement precision to the Heisenberg limit. But 
the measurement of the parity P is very difficult for sys- 
tems of large N because it requires the ability to distin- 
guish between odd and even numbers of spin-down par- 
ticles. Combining ideas from precision spectroscopy and 
quantum information processing, it has demonstrated the 
Heisenberg-limited interferometry in an experiment with 
three beryllium ions ^ISQi] , in which the nonlinear interac- 
tion from the laser-induced two-ion transitions has 
been used to prepare the NOON state and transfer the 
phase information into amplitude information. However, 
these protocols seems very difficult to scale to many-body 
systems of large N . 

For the systems of Bose condensed atoms, there are 
several methods for preparing high NOON states. Sim- 
ilar to prepare the spin squeezing, the high NOON 
states c an be per iodically generated by in the time evo- 
lutions |l65l - ll67j |. For double- well BEC systems, the 
macroscopic superposition states of all atoms in one well 
(a type of NOON st ates) could be generated by phase 
engineering p8l - [i7nl| . The efficient creation of NOON 
states has also b een demonstrated by utilizing quantum 
tunnelling |l7ll Il72l |. Moreover, it has also suggested 
that the high NOON states can be generated by an ef- 
fective interaction between two atoms from coherent Ra- 
man proces ses t hrough a two-atom molecular interme- 
diate state [l73j |. Recently, it finds that the tunneling 
resonances in a titled double-well BEC system lead to 
a dynamical scheme for creating few -body NOON-like 
macroscopic superposition states |l74| . Below, we focus 
on discussing how to prepare high NOON states with adi- 
abatic evolutions and then use them to accomplish a MZ 
interferometry [11]. 

Based upon the quantized BJJ, which is described 
by the giant-spin Hamiltonian (61) with D < Q (i.e. 
Ec < 0), one of us suggest an alternative protocol for 
preparing high NOON states and then performing the 
Heisenberg-limited MZ interferometry via adiabatic evo- 
lutions through critical points [ll| . In principle, this pro- 
tocol is suitable for many-body systems of thousands of 
Bose condensed atoms. In analogy to an optical MZ in- 
terferometer, the two collective spin states j^yi —-y) and 
-f^) are used as the two paths, and the beam split- 
ters are realized by the adiabatic processes through the 
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bifurcation point, see Fig. 10. Below, we give some key 
points of this protocol. 



Input I 





2 2/ 
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N N\ 
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FIG. 10: Schematic diagram for the adiabatic MZ interferom- 
etry via a quantized BJJ described by the giant-spin Hamilto- 
nian (61). The input state is a SU(2) coherent spin state along 
the X-axis. The beam splitters BSl and BS2 are achieved by 
the adiabatic processes through the critical point for a bifur- 
cation between normal and self-trapped states. 

The initial state is a CSS along the x-axis, which is the 
ground state for the giant-spin Hamiltonian of Bz = 0, 
> and D = 0. The desired NOON state is an su- 
perposition of two degenerated ground states for the case 
oi Bz — 0, Bx = and D < 0. It has demonstrated that 
a Hopf bifurcation from normal to self-trapped states 
takes place when the ratio \D/Bx\ increases from be- 
low to above the critical value A/N. This bifurcation, 
which relates to spontaneous symmetry breaking, can be 
looked as a continuous quantum phase transition. Keep- 
ing Bz = and D < unchanged, starting from the 
CSS along the x-axis, the NOON state can be adiabati- 
cally prepared by decreasing B^ from B^ ^ \D\ to zero, 
which accompanies with the occurrence of a dynamical 
bifurcation. This is to say, the first beam splitter could be 
achieved by the adiabatic evolution from \Bx/D\ 3> 1 to 
\Bx/D\ = 0. Even in the presence of Bz originated from 
the asymmetry 5, the NOON state can be still prepared 
with high fidelity ii S < Sc- 

After accomplishing the first beam splitter, the Joseph- 
son coupling between two modes is completely switched 
off, the state vector processes about the z-axis in 
the free evolution described by the unitary operator 

exp —iluqISz and accumulates a phase = cogt. Here, 

ujQ is the transition frequency between states \l) and |t)- 
Thus the state at the end of the free evolution is 



I*) 



exp [~i(t)Sz) ■ ( 



N_ 

2 ' " 



) + 



N_ 

2 ' ■ 



)) 



= ^ e+'f^|f,-f)+e-f^|f,+f) .(91) 



The second beam splitter is achieved by the inverse 
process of the first beam splitter. That is, keeping Bz ^ 
and D < unchanged, B^ is adiabatically increased from 
zero to Bx ^ \D\. After completing the second beam 
splitter, due to the dynamical symmetry breaking in the 
beam splitter, the final state is a superposition of the 
ground and first-excited states for the strong coupling 
limit Bx ^ II?! and the phase information of the state 
(91) are transferred into the amplitude information of the 



two lowest eigenstates, that is 



vE-)^ = cos ( -0 ) |GS) + e^" sin ( -<t> ] |FES) . (92) 



N 



Here, |GS) and |FES) denote the ground and first- 
excited states, respectively. In the strong coupling limit, 
\Bx/D\ > 1, the ground state |GS) = |f ,+f ) is a 
state of all spins along the positive x-axis direction, and 
the first-excited state |FES) — \^,+^- — l) is a state of 
{N — 1) spins along the positive x-axis direction and one 
spin along the negative x-axis direction. 

To estimate the phase 0, we measure the operator 

i?=|FES)(FES|-|GS)(GS|. (93) 

This leads to a measured signal and variance given by 



R)^ = - cos i^<t>) ' (^R) ^ = sin^ iN(f>) . 
Thus, the phase sensitivity is determined by 



A(f, 



d{R)^/d(j) 



1 



(94) 



(95) 



which reaches the Heisenberg limit 1/A^. 

The measurement of R requires the ability to distin- 
guish |GS) and |FES). However, in the strong cou- 
pling limit, it is difficult to directly distinguish these 
two eigenstates. So that one has to transfer these two 
states into distinguishable states. Fortunately, one can 
use the similarity of the two lowest eigenstates for sym- 
metric and asymmetric BJJs in the strong coupling limit 
and the non-degeneracy of the two lowest eigenstates for 
an asymmetric BJJ in the strong nonlinear limit. Af- 
ter accomplishing the second beam splitter, the system 
stays in the strong coupling limit, a proper asymmetry 
S is suddenly applied and then slowly decrease Bx from 
Bx ^ l-D] to zero. To keep the populations in |GS) and 
|FES) unchanged, the process must avoid the dynami- 
cal bifurcation and so that the asymmetry must satisfy 
\S\ < \Ea/2\ (the best choice is \S\ = \Ec/i\). Lastly, 
when Bx decreases to zero, |GS) = \N/2, —N/2) and 
|FES) = I A/2, +N/2), which can be easily distinguished. 

Comparing with the conventional Ramsey interferom- 
etry, in which the beam splitters are achieved by 7r/2 
pulses, in the interferometry protocol discussed above, 
the beam splitters are achieved by the adiabatic evolu- 
tions through the critical point for a Hopf bifurcation 
(also the critical point for a quantum phase transition). 
Moreover, in the adiabatic interferometry protocol, the 
beam splitters have important functions in addition to 
splitting and recombination. Besides splitting the in- 
put state, the first beam splitter prepares the desired 
NOON state for the Heisenberg-limited interferometry. 
Besides recombining the two paths, the second beam 
splitter transfers the phase information into amplitude 
information. 
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IV. CONCLUSIONS AND PERSPECTIVES 

In summary, we have given a review of recent progresses 
in nonlinear quantum interferometry via Bose condensed 
atoms. The nonlinearity originates from the atom-atom 
interaction, which can be tuned by applying external 
magnetic fields or lasers. Besides various exotic phe- 
nomena caused by the atom-atom interaction, the atom- 
atom interaction can be used to generate multi-particle 
entangled states, such as spin squeezed states and NOON 
states, which have been suggested for implementing high- 
precision measurements beyond the SQL provided by the 
non-entangled states. 

For the matter-wave interferometry based on the MQC 
of atomic BECs, the system can be well described by the 
MF theory. In ID systems, matter- wave solitons may be 
gradually generated in the interference process and more 
solitons will be generated for systems of more strong non- 
linearity. Similar to the SJ Js, the B JJs of linearly coupled 
BECs are widely used to detect and exploit the MQC of 
BECs. In addition to the Rabi oscillation, a direct signa- 
ture of MQC, the MQST and chaos may appear in BJJs 
of strong nonlinearity. By tuning the effective nonlinear- 
ity, the BJJ undergoes a symmetry-breaking transition 
from normal to self-trapped states, in which the criti- 
cal point depends on both the effective nonlinearity and 
relative phase. The slow passage dynamics across the 
critical point shows universal behaviors obeying the KZ 
mechanism for a continuous quantum phase transition. 

For the quantum interferometry based upon the many- 
body quantum coherence of Bose condensed atoms, the 
system must be described by the second quantization the- 
ory. Under the second quantization, a BJJ is described 
by two-mode Bose-Hubbard model, which is equivalent to 
an anisotropic giant-spin within external magnetic fields. 
The atom-atom interaction in this system is just the non- 
linear term of the one-axis twisting Hamiltonian for gen- 
erating spin squeezing. So that the atom-atom inter- 
action has been suggested for preparing spin squeezed 
states and then one can use those states for implement- 
ing high-precision interferometry beyond the SQL. More- 
over, by using adiabatic processes through the critical 
point, the NOON state can be prepared and then one 
can use this maximally entangled state for implementing 



Heisenberg-limited interferometry. It has suggested an 
adiabatic MZ interferometry, in which, the beam split- 
ters are realized by the adiabatic processes through the 
critical point and the two parts of the NOON states are 
used as two paths. 

It has been widely accepted that the multi-particle en- 
tanglement can enhance the measurement precision and 
the inter-particle interaction can be used to generate 
quantum entanglement. However, to design and build 
a practical quantum device based upon the quantum in- 
terferometry with multi-particle entangled states, there 
are still several important problems waiting to be solved. 

• Firstly, although only internal hyperfine or external 
motional states are used for interferometry, for an 
example an atom clock only use two internal states 
for interferometry. However, the atomic systems 
have both internal and external degre es of freedom 
and they may couple with each other [l75j | . There- 
fore, for interferometry with internal states, it is 
important to explore the effects from the mechani- 
cal motion of the center of mass. 

• Secondly, most of present works assume the sys- 
tems are closed systems. Actually, almost all prac- 
tical systems are open systems. Beyond the theo- 
retical frame for closed systems, we have to treat 
the systems as open systems and explore whether 
an interferometry protocol is robust against the en- 
vironment effects. Thus, it is vital to expl ore the 
effects from deco herence, temperature |l52j. dissi- 
pation [l76l Il77l| and particle losses [l78j etc. 

• Thirdly, all interferometry protocols involve dy- 
namical processes. It is essential to know whether 
imperfect excitations appear in the dynamical pro- 
cesses and how the imperfect excitations affect the 
interferometry. 
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